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ABSTRACT 

The  LINAC  at  NPGS ,  Monterey  was  used  to  accelerate  electrons 
to  energies  ranging  from  50  -  100  MeV.   These  were  used  to  study 

energy  losses  of  high  energy  electrons  in  aluminum  and  copper.   The 

2 
densities  of  each  material  ranged  from  0.7  to  2.8  gm/cm  . 

The  results  agreed  with  the  theory  of  Blunck  and  Westphal, 

unlike  previous  measurements  made  by  Breuer  who  found  disagreement 

between  his  experimental  results  and  his  interpretation  of  the  theory 

of  Blunck  and  Westphal,  particularly  at  energies  above  50  MeV  and 

with  aluminum  of  thickness  greater  than  that  yielding  a  density  of 

1  gm/cm  . 


UBRAF  . 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CALIF...  93940 


TABLE  OF  CONTENTS 

Section  Page 

I    INTRODUCTION  9 

II   THEORY  12 

III   EXPERIMENTAL  PROCEDURE  17 

IV   TREATMENT  OF  DATA  22 

V   RESULTS  31 

VI   DISCUSSION  61 

BIBLIOGRAPHY  63 

APPENDICES 

A.  Derivation  of  Expression  for  Recoil  Energy  65 

B.  Computer  Program  for  Counting  Rate  Corrections  67 
(Miller's  method) 

C.  Computer  Program  for  Counting  Rate  Corrections  68 
(rate  less  than  6  counts  per  second) 

D.  Computer  Program  for  Computing  Various  Distribution    69 
Parameters  (Aluminum) 

E.  Computer  Program  for  Computing  Various  Distribution 
Parameters  (Copper)  70 

F.  Computer  Program  for  Distribution  Calculations  72 


LIST  OF  TABLES 

Table  Title  Page 

I  DISTRIBUTION  CONSTANTS  13 

II  STERNHEIMER  CONSTANTS  15 

III  DISTRIBUTION  PARAMETERS  -  ALUMINUM  28 

IV  DISTRIBUTION  PARAMETERS  -  COPPER  29 

V  ENERGY  LOSS  IN  ALUMINUM  32 

VI  ENERGY  LOSS  IN  COPPER  33 


LIST  OF  ILLUSTRATIONS 

Figure           Title  Page 

1  Experimental  Arrangement  (top  view)  18 

2  Experimental  Arrangement  (side  view)  20 

3  Energy  Distribution,  74.63  MeV  Incident  23 
Energy,  Copper  Absorbers 

4  Histogram  Process  25 

5-16   Predicted  Curves  and  Observed  Points  35-46 
(Aluminum) 

17-28   Predicted  Curves  and  Observed  Points  47-58 
(Copper) 

29  Predicted  and  Observed  Half  Widths  (Aluminum)    59 

30  Predicted  and  Observed  Half  Widths  (Copper)  60 


ACKNOWLEDGMENTS 

The  author  wishes  to  express  his  appreciation  to  Professor 
Fred  R.  Buskirk  for  his  inspiration  and  outstanding  assistance 
in  carrying  out  this  project,  to  Professor  John  N.  Dyer  for  his 
invaluable  help  with  both  the  theoretical  aspects  and  the  actual 
performance  of  the  experiment,  and  to  Professor  Franz  A.  Bumiller 
who  initially  suggested  this  as  a  fruitful  area  of  investigation 
and  provided  overall  supervision. 

The  assistance  of  Mr.  F.  Ryan  and  2nd  Lt.  L.  P.  Gaby  in 
debugging  computer  programs  is  greatly  appreciated.   Valuable 
help  was  also  given  by  Lt.  J.  W.  Armstrong,  Lt .  N.  M.  Ferriter, 
and  Lt.  J.  W.  Stewart;  members  of  the  LINAC  group. 


I.   INTRODUCTION 

Electrons  passing  through  media  lose  energy.   The  principal 
losses  are  through  bremsstrahlung  (braking  radiation) ,  and 
ionization  and  excitation  of  the  atomic  electrons.   Due  to  the 
great  difference  in  the  masses  of  the  incident  electron  and  the 
nuclei,  virtually  no  energy  is  lost  to  recoil  of  the  nuclei. 

If  a  monoenergetic  beam  of  electrons  passes  through  a  slab 
of  matter,  the  transmitted  electrons  will  emerge  with  energy  less 
than  the  incident  beam.   In  addition  the  transmitted  electrons 
will  not  all  lose  the  same  amount  of  energy,  so  that  a  distribution 
of  energies  results.   If  a  polyenergetic  beam  passes  through  the 
slab,  it  will  already  have  an  energy  distribution  when  incident 
on  the  slab.   By  interaction  with  the  matter,  the  mean  of  the 
distribution  will  be  at  a  lower  energy  and  the  distribution  will 
be  broadened. 

From  theoretical  considerations  of  the  passage  of  electrons 
through  relatively  thin  materials  one  derives  a  most  probable  energy 
loss,  Q  ,  and  the  width  of  the  energy  distribution  of  the  exiting 
electrons.   The  width  considered  is  the  full  width  of  the  distri- 
bution at  half  the  maximum  amplitude.   This  will  be  called  the 
"half  width". 

When  the  incident  electrons  scatter  from  the  nuclei,  it  lengthens 
their  path  through  the  medium  (as  well  as  causing  bremsstrahlung) . 
For  lengthening  is  shown  to  be  insignificant  for  the  energies 
thicknesses  and  media  which  were  examined  by  Miller  (1) ,  who  based 
his  computations  on  the  work  of  Yang  (2)  with  data  from  Rossi  (3) . 
He  shows  that  for  incident  energies  greater  than  50  MeV  and  absorber 


thicknesses  less  than  that  yielding  5  grams  of  the  material 
per  square  centimeter,  the  predicted  path  increases  are: 
average  path  lengthening,  less  than  0.3%;  most  probable  path 
lengthening,  less  than  0.2%;  and  maximum  path  lengthening  for 
90%  of  the  electrons,  less  than  0.6%.   The  experimental  values 

to  be  reported  had  incident  energies  greater  than  53  MeV,  and 

2 
absorber  thicknesses  which  yielded  a  maximum  of  2.9  grams/cm  . 

Several  theoretical  treatments  of  these  energy  losses  (A, 
5,  6,  7,  8,  and  9)  and  a  few  measurements  (1,  10,  11,  12,  and 
13)  are  extant.   Breuer  (13)  found  agreement  with  the  theory 
of  Sternheimer  (9)  concerning  most  probable  energy  loss.   He 
found  half  widths  systematically  larger  than  the  values  calculated 
by  Blunck  and  Westphal  (4) . 

Miller  suggested  that  Breuer  erred  in  his  interpretation  of 
Sternheimer  when  he  concluded  that  the  most  probable  loss  due  to 

ionization  and  excitation,  Q  „,  varies  linearly  with  absorber 

pi  J 

thickness.   Miller  also  wished  to  examine  the  discrepancy  between 
the  half  widths  predicted  by  Blunck  and  Westphal  (B  &  W)  and  those 
observed  by  Breuer. 

Consequently  Miller  used  the  linear  electron  accelerator  (LINAC) 
at  the  Naval  Postgraduate  School  to  measure  the  energy  loss  of  electrons 
passing  through  aluminum,  with  incident  energies  between  60  and  100 
MeV.   Unfortunately  many  of  his  measurements  were  inconclusive  due 
to  geometrical  problems  encountered  in  the  gathering  of  data. 

Therefore  these  measurements  were  made  again  with  a  different 
experimental  arrangement.   Then  a  similar  set  of  data  was  taken  using 
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copper  as  the  absorbing  medium.   These  data  (for  both  aluminum 
and  copper)  agree  with  the  energy  loss  predictions  of  Blunck 
and  Westphal. 
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II.   THEORY 

Blunck  and  Westphal  (B  &  W)  (A)  have  developed  a  theoretical 

energy  distribution  for  the  exiting  electrons  which  considers 

losses  due  to  ionization,  excitation  and  radiation.   It  is  based 

on  the  works  of  Landau  (5) ,  Blunck  and  Leisegang  (6) ,  Eyges  (7) , 

and  Bethe  and  Heitler  (8) .   Their  distribution  is  expressed  in 

terms  of  a  dimensionless  parameter  A,  related  to  the  energy  loss 

Q  by: 

n   -r-    In  E. 

*- ^ -ar^  -  Li" 

Q   is  the  average  energy  loss  through  ionization  (in  MeV) . 

E.  is  the  incident  energy  of  the  electrons  (in  MeV).  and  R  is  the 
l 

thickness  (in  cm)  of  the  absorbing  medium.   The  constant  a  is  a 

function  of  the  atomic  number  Z,  the  atomic  weight  A,  and  the 

density  in  grams  per  square  centimeter  p  of  the  absorbing  medium; 

and  the  velocity  (6  =  v/c)  of  the  electron.   One  finds  from  B  &  W 

that: 

n    iq/    z           P                                fKev\ 
a  =  0.154  —     x  — r-  (  

A  g2  \   cm/ 

The  probability   distribution  W(Q)dQ   is   given  as: 

/     \  aR 
W(Q)dQ   =    B   x    y^j         x   FaRjb2       (A)    dA 


and 


a  -  An) 


F^v  k2    <a)    H   T(aR+  l)(y)  Z       /     2  l^F D-ar      >-o ^ 

«R,b  2  n     ^  +  Yn)l        aR  \/l2  +  Y   2 
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B  is  a  normalizing  constant.   The  quantities  G  .  Y  ,  X 

n   n   n 

are  constants  given  for  appropriate  values  of  the  summation 
constant,  n,  in  Table  I  below. 

TABLE  I.   DISTRIBUTION  CONSTANTS 


n  1 


C  0.174 

n 

Y  0.0 

n 

X  1.8 

n 


a  is  a  constant  related  to  the  losses  due  to  radiation.   It 
was  taken  by  B  &  W  from  Bethe  and  Heitler  and  is  given  as: 


2 

a  =  1.40  x  10"3  x  f-  x  p[4/3  In  (183/  1/3)  +  1/g]    (— ) 

A  Z  cm 

2 
The  term  b   is  from  Blunck  and  Leisegang  and  is  given  by: 

9     9  IN  A2E. 

b  =   F-)  (1.5)   Z   S„S   x  In' 


2 

3 

4 

0.058 

0.019 

0.007 

3.0 

6.5 

11.0 

2.0 

3.0 

5.0 

s    Z        VlsU  -  B2)) 

where  the  summation  is  over  the  s  different  ionization  potentials  of 

the  atomic  electrons  of  the  absorber.   N  is  the  number  of  electrons 

s 

per  atom  or  molecule  of  the  absorber  with  the  ionization  potential  I  . 

s 

2 
(When  calculating  b   for  the  copper  absorbers,  the  ionization  potentials 

given  by  Sternheimer  (9(a))  were  used.   In  this  case,  the  factor  1.5 

becomes  1.3.   Another  useful  source  of  these  potentials  is  Bearden 

and  Burr  (14) . 


is  the  parabolic  cylinder  function  of  the 


indicated  arguments. 
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B  &  W  calculated  values  of  the  function  F  _  .  2(a)  for  4  values 

aR,b 

2 
of  b  ;  0,  3,  6,  and  9.   These  are  presented  graphically.   Curves 

were  drawn  for  aR  values  of  0,  0.05,  0.10,  0.15,  0.20,  and  0.25 

2 
for  each  value  of  b  .   They  were  plotted  as  functions  of  A,  with  A 

ranging  from  approximately  -7  to  +15.   For  a  given  incident  energy 

and  absorber  thickness,  the  distribution  W(Q)  is  directly  proportional 

to  F   ,2(A).   All  other  terms  are  normalization  terms.   Because  Q 

is  linearly  proportional  to  A,  one  can  write: 

Q  =  (aR)A  +  Q^ 
Q  "  Q. 


A 


aR 


where  Q   is  the  energy  loss  corresponding  to  A  equals  zero,  and  is 
o 

given  by: 

Qo  =  Q  -  aR  x  In  f  ^  )  +  1.116(aR) 

In  evaluating  Q,  the  average  energy  loss  due  to  ionization  and 
excitation,  the  equations  of  Sternheimer  (9)  are  used.   These  consider 
the  density  effect  due  to  polarization  of  the  medium.   His  equations 
for  the  average  energy  loss  can  be  combined  and  written  as: 


A  t 

Q-  A- 


?  ms 

B  +  0.43  +  In  E.  -  B  -  C  -  a  (Xn  -  log1A  £-) 
s  i  sv  1     &10  mc 


The  terms  A  ,  B  ,  C,  X. ,  a  ,  and  m  are  constants  characteristic 
s   s      1   s       s 

of  the  absorbing  material.   They  are  found  in  Sternheimer  9(c)  and 
are  tabulated  below  for  aluminum  and  copper. 
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TABLE 

II 

A 
s 

B 
s 

C 

A£ 

0.0740 

16.77 

-4.21 

Cu 

0.0701 

15.09 

-4.74 

X 


m 


Is        s 
3     0.0906    3.51 


0.119 


3.38 


From  the  relativistic  mass/energy  relationship  one  can  find  that 


P_  =  ££  =\ 

mc     ? 

mc 


2 T~ 
E  +  2mc  E 


mc 


Sternheimer  (9,b)  also  gives  an  expression  for  the  most  probable 

energy  loss  due  to  ionization  and  excitation,  Q  . .   Of  interest  is 

pi 

a  comparison  of  Q  . ,  the  most  probable  loss  from  B  &  W,  and  experimental 
results.   This  will  show  the  relative  effects  of  ionization  and  excit- 
ation compared  to  that  of  radiation  losses  on  the  most  probable  energy 

loss.   The  expression  for  Q  .  is: 

pi 


A  t 
s 


Pi 


B  +  1.06  +  In 
s 


s 


m 

-  C  -  a   (X.  -  login  ^-X 
s  V  1      10  mc 


If  one  assumes  a  given  E.  and  small  energy  losses,  this  expression 
reduces  to : 


Q  ,  -  K  t  +  K  t(ln  t) 


where  K  and  K_  are  constants  determined  by  the  previous  expression. 

This  shows  the  variation  of  the  most  probable  energy  loss  (due 

to  excitation  and  ionization)  with  absorber  thickness.   It  is  clearly 

non-linear  in  t,  contrary  to  the  interpretation  of  Breuer.   Dividing 

through  by  t  yields: 

Q   /t  =  K.  +  K_  (In  t) 
pi      1    2 
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This  shows  that  given  an  initial  energy  the  quantity  Q  ,/t  is  not 

uniquely  determined,  but  rather  varies  with  the  thickness  of  the 

absorber,  t.   Plotting  Q  . /t  versus  E,  results  in  a  family  of  curves, 

each  for  a  different  absorber  thickness,  as  shown  by  Miller  in  his 

figure  1;  this  is  not  a  single  straight  line  as  showrt  by  Breuer  in 

his  figure  4. 

Miller  shows  in  figure  3  a  series  of  curves  from  which  one  can 

derive  B  &  W's  prediction  for  the  values  of  most  probable  energy  loss 

and  half  width.   He  also  plots  Q   versus  t  for  an  E.  of  70  MeV  as 

pi  l 

figure  2.   Figure  3  shows  Q  ,  the  most  probable  energy  loss  when 
radiation  losses  are  also  considered,  versus  t  for  70  MeV. 
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III.   EXPERIMENTAL  PROCEDURE 

In  order  to  measure  the  energy  loss  of  electrons  passing 
through  matter,  the  LINAC  was  used  as  a  source  of  the  high  energy 
electrons.   The  energy  of  the  exiting  electrons  was  measured  by 
using  the  16  inch  magnetic  spectrometer  described  by  Kenaston,  Luke, 
and  Sones  (15) .   The  energy  of  the  incident  beam  was  measured  by 
using  a  nuclear  magnetic  resonance  probe  to  measure  the  field  of  the 
deflection  magnet  of  the  LINAC. 

If  the  spectrometer  had  been  placed  in  a  zero  degree  position, 
any  beam  current  large  enough  to  be  stable  would  have  saturated  the 
counting  system.   For  this  reason,  a  3.3  mil  aluminum  foil  was  placed 
at  an  angle  of  45   to  the  electron  beam.   This  was  used  to  scatter  out 
from  the  main  beam  a  secondary  flux  of  electrons.   Close  behind  the 
scattering  foil  were  the  absorbers  (see  figure  1) .   These  were  parallel 
to  the  beam  line.   The  spectrometer  was  set  at  an  angle  of  90   to  the 
beam  approximately  30  cm  from  the  absorbers. 

Miller's  initial  set-up  was  the  same  except  that  his  absorbers 
were  22  cm  from  the  scattering  foil,  mounted  on  a  rotating  wheel.   In 
his  "Discussion"  he  shows  that  this  lead  to  a  broadening  of  the  energy 
distributions  of  the  exiting  electrons.   This  was  due  to  electrons 
which  were  scattered  by  the  absorber  into  the  spectrometer  which 
resulted  in  their  appearing  to  have  come  from  a  displaced  source  with 
different  energy.   This  is  discussed  at  length  by  Miller.   Another 
slight  difference  is  that  he  set  the  spectrometer  at  45  . 

For  this  reason  the  absorbers  were  mounted  close  to  the  scattering 
foil  (7  mm  away  at  one  end  and  39  mm  at  the  other.)   In  order  to  position 
absorbers  of  different  thicknesses  in  place,  a  laminate  of  differing 
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lengths  of  absorber  material  was  raised  and  lowered  by  a  piston. 
Figure  2  illustrates  this  arrangement.   Because  the  target  chamber 
was  evacuated  to  approximately  50  microns  there  was  no  air  between 
the  layers  of  absorber  material  to  affect  the  results.   The  absorbers 
were  parallel  to  the  beam  and  2.2  cm  away  from  the  beam  center  line. 
They  did  not  scatter  electrons  from  the  main  beam  into  the  spectrometer. 
The  slot  over  which  no  scattering  foil  was  placed  was  used  to  measure 
the  experimental  background.   The  topmost  position  which  had  scattering 
foil  but  no  absorber  behind  was  used  to  measure  the  energy  distribution 
of  the  electrons  incident  on  the  absorbers. 

The  vacuum  coupling  described  by  Miller  was  used  to  couple  the 
spectrometer  to  the  target  chamber,  eliminating  all  windows  between 
the  absorbers  and  the  spectrometer  magnet. 

The  spectrometer  and  counting  system  measure  the  number  of  electrons 
within  a  small  energy  interval  entering  the  magnet.   For  this  reason 
many  measurements  at  slightly  different  energies  are  necessary  to 
construct  the  profile  of  the  energy  distribution.   In  order  to  ensure 
that  each  measurement  represented  the  same  amount  of  secondary  beam 
incident  on  the  absorbers  a  secondary  emission  monitor  (SEM)  was  used 
to  monitor  the  main  electron  beam.   The  output  of  the  SEM  was  measured 
by  a  Cary  integrating  electrometer,  with  a  10  microfarad  input  capacitor. 
When  this  reached  a  predetermined  charge  counting  was  stopped.   In  this 
manner  the  data  were  standardized  to  a  given  amount  of  charge  from  the 
beam  incident  on  the  scattering  foil.   Typical  values  used  are  100 
millivolts  to  1  volt  which  represent  10  to  100  microcoulombs  of  beam 
charge  (because  of  the  10%  efficiency  of  the  SEM) . 
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FIGURE  2.   Experimental  Arrangement  (side  view) 
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Data  taken  were:   spectrometer  energy  setting,  number  of  counts 
observed,  time  of  integration,  and  voltage  to  which  the  charge  on 
the  capacitor  was  integrated.   At  approximately  1  MeV  intervals  the 
background  was  measured  by  integrating  charge  with  the  closest  (to 
the  beam)  absorber  in  place  but  no  scattering  foil.   These  values 
were  used  later  in  determining  how  much  background  to  subtract  from 
the  number  of  counts  observed. 
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IV.   TREATMENT  OF  DATA 

The  theory  presented  earlier  presumes  a  monoenergetic  source 
of  electrons.   While  the  LINAC  produces  an  essentially  monoenergetic 
beam  of  electrons  there  is  a  finite  width  distribution  in  energy 
of  the  electrons  comprising  the  beam.   The  secondary  beam,  which  is 
incident  on  the  absorbers,  is  scattered  from  the  main  beam  by  a  thin 
aluminum  foil  and  this  scattering  adds  to  the  width  of  the  energy 
distribution.   Thus  the  secondary  beam  incident  on  the  absorbers  had 
an  experimentally  measured  average  half  width  of  0.5%  of  the  incident 
energy.   This  width  cannot  be  considered  monoenergetic.   A  plot  of 
the  distribution  in  energy  of  the  electrons  incident  on  the  absorbers 
and  the  resulting  distributions  in  energy  of  the  electrons  after 
passing  through  the  several  absorbers  is  shown  as  figure  3.   This 
is  a  plot  of  the  values  observed  with  incident  energy  74.63  MeV  and 
copper  absorbers.   It  is  typical  of  all  the  other  combinations  of 
energy  and  absorber  material. 

The  difficulty  in  deriving  a  predicted  energy  loss  and  half  width 
stems  from  the  fact  that  neither  is  the  incident  beam  monoenergetic,  nor 
is  there  a  simple  mathematical  expression  for  its  distribution.   If  it 
were  monoenergetic,  the  resulting  distributions  after  passing  through 
the  absorbers  could  be  derived  from  the  theory  presented  in  section 
II  in  a  straightforward  manner.   If  the  beam  were  polyenergetic  but 
were  distributed  in  a  simple  manner,  the  resulting  distributions  could 
also  be  derived  from  the  theory  in  a  not  so  straightforward  but  soluble 
manner.   One  would  perform  the  integration: 


R(Ef)  =   J  W(Ei  -  Ef)  x  D(Ei)  dEi 

E, 

l 
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Figure  3.  Energy  Distribution,  74.63  MeV  Incident  Energy  Copper  Absorber* 
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where  R(Ef)  is  the  resulting  distribution,  W(Ei  -  Ef)  is  from  the 
theory,  and  D(Ei)  is  the  hypothesized  known  incident  distribution 
function.  The  range  of  integration  would  be  over  the  range  of 
incident  energies.   This  would  be  a  difficult  integration,  at  best. 

A  very  good  approximation  is  to  use  the  histogram  method 
described  by  Miller.   The  incident  distribution  is  divided  into  a 
large  number  of  equal  width  energy  slices  and  each  of  these  is 
considered  a  monoenergetic  source  with  amplitude  given  by  the  value 
of  real  distribution  at  the  mid-point  of  the  slice.  The  theoretical 
distribution  was  similarly  divided  into  slices  with  the  same  width  in 
energy  as  the  incident  beam  distribution  slices. 

Then  the  amplitudes  of  the  observed  incident  distribution  and 
theory  curves  from  B  &  W  were  fed  into  the  IBM  360  computer  at  the 
Naval  Postgraduate  School  and  handled  in  the  following  manner.   Each 
experimental  amplitude  was  multiplied  in  turn  by  each  of  the  amplitudes 
from  the  theory  curves  and  the  resulting  products  were  stored  in  bins 
together  with  the  appropriate  energy  value  for  that  bin,  determined  by 
the  theory  of  B  &  W.   Succeeding  products  were  added  to  the  contents 
of  the  appropriate  bins.   In  this  manner  resulting  energy  distribution 
predictions  were  generated.   This  process  is  illustrated  in  figure  4. 
The  observed  values  for  each  absorber  and  energy  combination  had 
previously  been  plotted,  and  a  smooth  curve  drawn  through  the  points. 
The  predicted  curves  were  normalized  by  the  computer  to  have  the  same 
maximum  amplitude  as  the  faired  in  observed  value  curves. 

Figure  4  shows  two  histograms  D(E)  and  W(Q)  being  used  to  generate 
a  third,  R(E).   They  correspond  to  the  observed  initial  distribution, 
B  &  W  curves,  and  predicted  distribution,  respectively. 
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2 
The  B  St   W  curves  used  were  those  for  b  =0.   Miller  shows 

(c.f.  his  figure  8)  that  there  is  no  significant  difference 

2  2 

between  the  curve  predicted  by  the  b  =  0  curves  in  B  &  W  and  b  =3. 

2 
In  this  work,  b   ranged  from  a  low  of  0.088,  corresponding  to  the 

2 
thickest  aluminum  absorber  (2.86  gms/cm  )  and  53.6  MeV;  to  a  high  of 

2 
0.934,  corresponding  to  the  thinnest  copper  absorber  (0,71  gms/cm  ) 

2 
at  94.3  MeV.   From  this  one  can  see  that  the  b  =  0  curves  were 

appropriate. 

The  B  &  W  curves  are  truncated  at  A  =  15.   This  results  in 

minor  discrepancies  at  the  low  energy  tail  of  the  predicted  distribution. 

Therefore  they  were  extended  exponentially  to  A  =  30;  by  matching  the 

amplitude  and  slope  at  A  =  15.   This  would  not  affect  the  predicted 

half-width  but  should  mojfe  accurately  represent  the  low  energy  tail. 

In  one  case  it  became  necessary  to  extrapolate  from  aR  =  .25  to  obtain 

-2 
the  curve  for  aR  =  .30;  for  the  Cu  target  with  t  =  2.85  gm-cm 

Incident  energies  used  in  calculating  other  parameters  were 

calculated  by  taking  the  value  measured  by  the  nuclear  magnetic 

resonance  probe  and  subtracting  from  it  energy  losses  due  to  recoil; 

and  ionization,  excitation,  and  radiation  when  passing  through  the  thin 

scattering  foil.   Recoil  losses  were  calculated  from  the  formula: 


E    / 2  sin  8/2 
T  =   °    ZE    T~ 

M  l     \  1  +  — |  sinZ6/2 
c   \    M  2 

c 

from  Hofstadter  (16),  and  checked  by  the  relationship  derived  in 
Appendix  A.   In  the  Hofstadter  formula  E  is  incident  electron  energy, 
M  the  mass  of  the  aluminum  nucleus,  and  9  was  90  in  this  case.   (This 
is  not  inconsistent  with  the  opening  remarks  that  in  passage  through 
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the  absorber,  the  losses  to  recoil  are  insignificant.   Very  few  of 
the  electrons  incident  on  a  slab  are  scattered  out  to  an  angle  of 
90°,  yet  virtually  all  of  the  electrons  incident  on  the  absorbers 
have  been  scattered  to  an  angle  near  90°  by  the  foil.   Because  of 
the  great  difference  in  the  masses  the  losses  are  small,  typically 
0.3  MeV  loss  at  94.5  MeV,  but  can't  be  neglected.)   The  values  for 
ionization,  excitation,  and  radiation  losses  were  taken  from  Berger 

and  Seltzer  (17) . 

The  B  &  W  theoretical  values  reported  in  the  Results  section 
are  not  the  Q  from  the  B  &  W  theory  because  the  B  &  W  curves  do 
not  peak  at  A  =  0  and  the  predicted  most  probable  energy  loss  is 
slightly  greater  than  Q  .   The  reported  values  were  obtained  by 
subtracting  the  energy  value  of  the  peak  of  the  predicted  curve 
from  that  of  the  peak  of  the  incident  beam.   The  predicted  half 
widths  were  measured  from  the  predicted  curves.   These  predicted 
curves  and  the  corresponding  experimental  points  are  shown  as  figures 
5  through  28  in  the  Results  section. 

Distribution  parameters  for  aluminum  and  copper  are  shown  below 

as  Tables  III  and  IV. 
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TABLE  III   DISTRIBUTION  PARAMETERS -ALUMINUM 

t(gm/cm2)      Ei(MeV)  b2          aR(MeV)          R(cm) 

0.730         96.93  0.374        0.0542          0.2743 

74.63  0.361        0.0542 

53.57  0.345        0.0542 

1.441         96.93  0.190        0.1070          0.5418 

74.63  0.183        0.1070 

53.57  0.175        0.1070 

2.146         96.93  0.127        0.1593          0.8067 

74.63  0.123        0.1593 

53.57  0.117        0.1593 

2.859         96.93  .  0.096        0.2123          1.0749 

74.63  0.092        0.2123 

53.57  0.088        0.2123 

Material:   EC  grade  aluminum,  purity  99.4%+ 
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TABLE  IV  DISTRIBUTION  PARAMETERS-COPPER 

2  2 

t(gm/cm  )      E  (MeV)        b  aR(MeV)  R(cm) 

0.711         94.30        0.934        0.050  0.080 

74.76        0.904        0.050 
52.84        0.858        0.050 

1.423         94.30        0.467        0.09999         0.160 
74.76        0.452        0.09999 
52.84        0.429        0.100 

2.134         94.30       0.311        0.14999         0.240 
74.76        0.301        0.14999 
52.84        0.286        0.150 

2.845         94.30        0.233        0.19998         0.320 
74.76        0.226        0.19999 
52.84        0.214        0.19999 

Material:   Hard  drawn  copper,  purity  99.0%+ 

The  values  presented  as  measured  or  observed  number  of  counts  are 
the  number  obtained  by  taking  the  raw  data,  correcting  it  for  the 
effect  of  losses  due  to  high  counting  rate,  and  subtracting  off  the 
background.   For  aluminum,  a  method  similar  to  that  used  by  Miller 
was  used  to  correct  for  high  counting  rate.   His  relationship  was: 

N  =  EN(1  +  (S)  Lfo 

where  N  is  the  corrected  number,  EN  is  the  raw  data,  and  T  is  the  time 

EN 
taken  to  count  EN  counts.   The  term  — ■=:   is  the  counting  rate  at  which  the 
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3.3 
number  EN  was  counted.   The  factor  TT7X  was  derived  by  varying  the 

counting  rate  with  a  standard  set  up  and  observing  the  response  of 

the  counters.   This  response  was  graphed  (counts  versus  counting 

rate)  and  the  slope  of  this  graph  determined  the  numerical  value  of 

the  correction  factor. 

The  relationship  used  for  aluminum  in  this  case  was: 

F      FN 
N  =  EN  J  (1  -  (  Yooo  X  T^  "  Back8round 

F 
Here   ..  ■„.  is  the  numerical  value  obtained  by  graphing  the  response 

of  the  counters  after  all  data  was  taken  at  each  energy. 

For  copper  the  beam  was  kept  down  to  a  value  which  would  produce 

a  maximum  of  6  counts  per  second.   When  this  is  done,  the  correction 

becomes : 

r  =  r  (1  +  .6r  ) 
to        o 

where  r  is  the  true  number  of  counts  per  machine  pulse,  and  r  is  the 
t  o 

number  of  counts  per  machine  pulse  observed.   This  requires  converting 
the  observed  counts  and  time  to  counts  per  machine  pulse,  but  as  shown 
by  Browman  (18)  it  will  yield  more  accurate  results  because  it  is  not 
dependent  on  the  instantaneous  shape  of  the  machine  pulses.   Computer 
programs  for  correcting  the  observed  counts  are  included  as  Appendices 
B  and  C. 
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V.   RESULTS 

The  values  predicted  by  theory  and  experimental  results  of  most 

probable  energy  loss  and  distribution  half  widths  for  both  aluminum 

and  copper  are  shown  in  tables  V  and  VI .   The  column  headed  S  under  Q 

P 

is  the  Sternheimer  prediction,  B  &  W  is  the  Blunck  &  Westphal  prediction, 

As  one  can  readily  see,  the  measured  data  fit  the  B  &  W  theory 
well,  and  differ  significantly  from  the  Sternheimer  predictions.   When 
considering  energy  loss  in  aluminum,  the  observed  data  differ  from  the 
B  &  W  curves  by  an  average  of  2  1/4%  and  differ  from  Sternheimer  by  7%. 
With  copper,  the  difference  is  more  dramatic;  average  difference  from 
B  &  W  is  1  1/2%,  average  difference  from  Sternheimer  is  8  1/3%.   The 
half -widths  observed  agree  with  the  B  &  W  theory  also.   Of  the  24 
half -widths  measured,  6  fall  further  away  from  the  predicted  value  than 
the  error  limits  assigned  to  the  observed  points.   Because  of  the  way 
the  half  width  predictions  were  determined,  error  limits  should  be 
assigned  to  the  predicted  values.   These  come  from  the  uncertainty  in 
the  exact  energy  of  the  incident  beam.   If  a  +  0.3%  value  is  used  as 
the  uncertainty  in  energy  (see  Oberdier  (19)  and  Miller  (1)),  then 
this  would  result  in  only  three  of  these  points  being  out  of  the  error 
limits.   This  is  quite  good  agreement.   If  there  is  a  systematic  error 
it  is  that  the  observed  half  widths  are  narrower  than  those  predicted 
and  is  the  opposite  effect  from  that  observed  by  Breuer. 

In  order  to  determine  the  tolerances  in  the  observed  energy  losses 
and  half  widths,  error  bars  were  put  on  the  observed  points,  and  extreme 
case  curves  were  plotted. 
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TABLE  V  ENERGY  LOSS  IN  ALUMINUM 


Q  (MeV) 

HW 

(MeV) 

E . (MeV) 
a. 

2 
t(gm/cm  ) 

B  &  W 

S 

Experiment 

Theory 

Experiment 

96.93 

0.730 

1.082 

0.977 

1.17 

0.67 

0.66 

1.441 

2.166 

2.002 

2.16 

1.00 

1.01 

2.146 

3.247 

3.044 

3.24 

1.32 

1.19 

2.859 

4.370 

4.117 

4.32 

1.79 

1.48 

74.63 

0.730 

1.066 

0.976 

1.04 

0.59 

0.51 

1.441 

2.094 

2.000 

2.18 

0.91 

0.91 

2.146 

3.174 

3.041 

3.17 

1.21 

1.08 

2.859 

4.297 

4.113 

4.26 

1.61 

1.38 

53.57 

0.730 

1.051 

0.974 

1.03 

0.47 

0.49 

1.441 

2.088 

1.996 

2.16 

0.75 

0.75 

2.146 

3.172 

3.036 

3.18 

0.99 

1.02 

2.859 

4.347 

4.106 

4.21 

1.35 

1.32 
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TABLE  VI   ENERGY  LOSS  IN  COPPER 


% 

(MeV) 

HW(MeV) 

E . (MeV) 

2 
t(gm/cm  ) 

B  &  W 

S 

Experiment 

Theory 

Experiment 

94.30 

0.711 

1.014 

0.840 

0.95 

0.69 

0.74 

1.423 

1.947 

1.750 

1.94 

1.12 

0.99 

2.134 

2.956 

2.685 

2.93 

1.62 

1.47 

2.845 

4.082 

3.637 

4.00 

2.31 

1.75 

74.76 

0.711 

0.913 

0.839 

0.92 

0.62 

0.66 

1.423 

1.894 

1.747 

1.91 

1.01 

1.07 

2.134 

2.903 

2.681 

2.92 

1.51 

1.43 

2.845 

3.978 

3.633 

3.98 

2.16 

1.98 

52.84 

0.711 

0.950 

0.837 

0.94 

0.50 

0.52 

1.423 

1.880 

1.743 

1.87 

0.82 

0.84 

2.134 

2.886 

2.675 

2.87 

1.27 

1.16 

2.845 

3.959 

3.624 

3.80 

1.91 

1.75 
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Plots  of  the  predicted  curves  and  observed  points  are  included 
as  figures  5  through  28.   These  were  plotted  by  the  IBM  360  computer, 
A  plot  of  predicted  and  observed  half  widths  with  the  aluminum 
absorbers  is  figure  29;  copper,  figure  30. 
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Figure  5.  Predicted  Distribution,  Al 
t=0.730  gm/cm2   Ei=  96.93  MeV 
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Figure  6.  Predicted  Distribution,  Al 
t=  1.441  gm/cm    Et=  96.93  MeV 
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Figure  7.  Predicted  Distribution,  Al 
t=  2.146  gm/cm2   E^  96.93  MeV 
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Figure  8.  Predicted  Distribution  Al 
t=  2.859  gm/cni    Ef  96.93  MeV 
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Figure  10.  Predicted  Distribution,  Al 
t=  1.441  gm/cm2   Ei=  74.63  MeV 
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Figure  11.  Predicted  Distribution,  Al 
t=  2.146  gm/cm2   Ei=  74.63  MeV 
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Figure  12.  Predicted  Distribution,  Al 
t=  2.859  gm/cm2   E±=   74.63  MeV 


42 


E  (MeV) 

Figure  13.  Predicted  Distribution,  Al 
t=  0.730  gm/cm2   Et=   53.57  MeV 
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Figure  14.  Predicted  Distribution,  Al 
t=  1.441  gm/cm2   Et=  53.57  MeV 
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Figure   15.    Predicted  Distribution,   Al 
t=  2.146  gm/cm2       E±=  53.57  MeV 
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Figure    16.    Predicted  Distribution,   Al 
t=  2.85  9   gm/cm2       E.=  53.57  MeV 
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Figure  17.  Predicted  Distribution,  Cu 
t=  0.711  gm/cm2   E,-  94.30  MeV 
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Figure  18.  Predicted  Distribution,  Cu 
t=  1.423  gm/cm2   Ej.-  94.30  MeV 
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Figure  19.  Predicted  Distribution,  Cu 
t=  2.134  gm/cm2   Ei=  94.30  MeV 
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Figure  20.  Predicted  Distribution  Cu 
t-  2.845  gm/cm2   E±=   94<3()  MeV 
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Figure  21.  Predicted  Distribution,  Cu 
t=  0.711  gm/cm2   Ei=  74.76  MeV 
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Figure  22.  Predicted  Distribution  Cu 
t=  1.423  gm/cm2   Ei=  74.76  MeV 
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VI.   DISCUSSION 


An  attempt  was  made  to  match  absorber  thicknesses  and  energies 

with  Miller's  work.   The  thinnest  absorber  here  was  to  have  been 

2 
the  same  as  Miller's  second  thinnest  (0.730  versus  0.711  gm/cm  ); 

the  next  the  same  as  Miller's  fourth  thickest  (1.441  versus  1.522), 

etc.   The  differences  can  be  explained  by  the  inaccuracies  in 

machining  the  absorbers,  and  the  fact  that  a  different  density  for 

3 
aluminum  was  used.   Miller  used  the  tabulated  value,  2.69  gms/cm  . 

The  supplier  of  the  aluminum  indicates  an  approximate  density  of 

3 
2.66  gm/cm  .   Weighing  a  sample  then  measuring  it  with  a  micrometer 

yielded  2.66,  as  did  weighing  another  sample  and  measuring  the  volume 

of  water  displaced  by  it.   Hence  the  value  of  the  density  of  aluminum 

3 
used  was  2.66  gm/cm  .   The  density  of  the  copper  samples  was  measured 

3 
by  the  last  two  methods  and  found  to  be  8.89  gm/cm  . 

Generally  speaking  the  results  agree  with  the  theory  of  Blunck 

&  Westphal,  indicating  that  both  the  theory  and  the  experimental 

technique  are  good.   They  also  agree  substantially  with  the  results 

obtained  by  Miller  for  energy  loss.   For  instance  Miller  found  at  53.7 

MeV  and  t  =  0.711  an  energy  loss  of  1.05  +  .07,  here  we  find  1.03  +  .03 

2 
at  53.57  MeV  and  t  =  0.730  gm/cm  .   Breuer  obtains  approximately  1  MeV 

2 
loss  at  53.6  MeV  with  absorber  thickness  =  0.73  gm/cm  .   (He  states 

-2 
AE  =1.36  MeV/gm  cm  )  but  indicates  that  with  higher  energies  and 
P 

thicker  absorbers  his  experimental  results  diverge  by  increasing 
amounts  from  the  theory.   The  conclusions  to  be  drawn  from  the  present 
results  tend  to  confirm  B  &  W  theory  and  indicate  that  Miller  is  correct 
in  his  assertion  that  Breuer  misinterpreted  it. 
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Variation  from  the  B  &  W  predictions  tend  to  be  larger  with 
thicker  absorbers.   Not  only  is  the  theory  weaker  in  these  regions 
(energy  losses  are  larger  compared  to  incident  energy) ,  but  also 
the  observed  values  are  most  inaccurate.   The  peaks  are  broad  and 
not  well  defined.   The  percentage  error  increases  as  the  number  of 
counts  decreases;  and  only  by  counting  for  a  prohibitively  long 
time  can  it  be  significantly  improved;  and  this  would  be  nullified 
to  some  extent  by  drift  in  the  operating  characteristics  of  the 
LINAC.   The  background  is  a  larger  fraction  of  the  total  number  of 
counts  observed,  and  errors  made  in  determining  the  proper  back- 
ground are  more  significant.   But  still  with  all  of  these  effects, 
the  experiment  and  theory  agree. 

Positioning  the  absorber  close  to  the  scattering  foil  has 
obviated  collimating  the  electrons  leaving  the  absorber.   This  was 
attempted  by  placing  a  lead  brick  with  a  hole  drilled  in  it  in  front 
of  the  spectrometer  entrance,  and  the  net  effect  was  an  unacceptable 
increase  in  the  background  with  no  commensurate  increase  in  energy 
resolution. 

Further  work  in  this  field  at  the  NPGS  LINAC  is  indicated,  using 
different  materials  as  absorbers.   One  would  expect  to  find  further 
confirmation  of  the  B  &  W  theory. 
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APPENDIX  A.   Derivation  of  Expression  for  Recoil  Energy.  (90° 
scattering) . 

E  =  Incident  electron  energy 

E  '  =  final  electron  energy 

M  =  rest  mass  of  the  struck  nucleus 

E  =  final  energy  of  the  struck  nucleus 
a 

From  conservation  of  mass  and  momentum: 

E  +  M  -  E  '  +  E   and  p2=p2+p'2 
e        e     a  e    *e 

from  these  we  find: 

E    '    =    (E     +  M  -  E   ) 
e  e  a 

E,2=(E     +M-E)2 
e  e  a 

E    ,2=    (E2+M2+E2+2EM-2EM-2EE) 
e  e  a  e  a  e  a 

2 
The  term  E  "  will  be  neglected  because  it  is  second  order  in  E  , 
a  a 

which  is  small  compared  to  E  and  M. 

E  ,2  =  E  2  +  M2  +  2E  M  -  2E  M  -  2E  E  . 
e      e  e      a      e  a 

2 
Similarly  p  is  small  compared  to  p   and  p   so  we  neglect  p  . 

,2     2 

P    =  P 
e      e 

2   v,2     2   _  ,2 

m  =  M  +  p   +  E 
e     e 

m?  =  M2  +  p  2  +  E  2  +  M2  +  2E  M  -  2E  (M  -  E  ) 
e     e  e      a      e 

from  which 

2  2  2  2 

M+p+E+M+2EM 

E      =  IS § ?_ 

a  2(M  +  E   ) 

e 


2  2  2 

2M     +  p       +  E       +  2E  M 

T     =  § § 5_     _  M 

a  2(M  +  E   ) 
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T      = 
a 


2M2  +  p   2  +  E  2  +  2E  M  -   2M2  -  2E  M 
^_e e e e 

2(M  +  E  ) 
e 


T     = 


2  JL  v  2  2  2  2  0      2    ,      2 

p+E  D+p+m  2p+m 

e  e  "e  e  e 


"a        2(M  +  E   ) 
e 


2(M  +  E   ) 
e 


2(M  +  E   ) 
e 


here  m  is  small  compared  to  the  other  terms. 

2 
Pe 
Ta  =  '(M  +  E  )       (wLthin  approximations  noted) 
e 
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APPENDIX  B.   Computer  Program  for  Counting  Rate  Corrections 
(Miller's  method) 
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Appendix  C.  Computer  Program  for  Counting  Rate  Corrections 
(rate  less  than  6  counts  per  second) . 
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APPENDIX  D.  Computer  Program  for  Computing  Various  Distribution 
Parameters  (Alumii;um 
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APPENDIX  E.  Computer  Program  for  Computing  Various  Distribution 
Parameters  (Copper). 
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APPENDIX  F.  Computer  Program  for  Distribution  Calculations 
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